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Thermoelectric effects in a double quantum dot system coupled to external magnetic/nonmagnetic 
leads are investigated theoretically. The basic thermoelectric transport characteristics, like ther- 
mopower, electronic contribution to heat conductance, and the corresponding figure of merit, have 
been calculated in terms of the linear response theory and Green function formalism in the Hartree- 
Fock approximation for Coulomb interactions. An enhancement of the thermal efficiency (figure of 
merit ZT) due to Coulomb blockade has been found. The magnitude of ZT is further considerably 
enhanced by quantum interference effects. Both the Coulomb correlations and interference effects 
lead to strong violation of the Wiedemann-Franz law. The influence of spin-dependent transport 
and spin bias on the thermoelectric effects (especially on Seebeck and spin Seebeck effects) is also 
analyzed. 
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I. INTRODUCTION 

Energy conversion based on thermoelectric properties 
of solid-state materials has been attracting recently much 
attention, especially in the case of nanoscale systems. Al- 
though thermoelectric phenomena have been known for 
long time, their efficiency in bulk materials was relatively 
low, and therefore systems of high thermoelectric effi- 
ciency are greatly desired. The efficiency is usually mea- 
sured by the dimensionless thermoelectric figure of merit 
ZT, ZT = ctS 2 T/k. Here, T stands for the temperature, 
a is the charge conductivity, S is the Seebeck coefficient 
(thermopower), and n is the thermal conductivity which 
generally includes the phonon and electronic contribu- 
tions. From this formula it is clear that to enhance the 
efficiency of a thermoelectric device one has to somehow 
enlarge the thermopower S and electric conductivity a, 
and reduce the thermal conductivity k. However, this 
task is difficult to be realized in conventional materials 
which obey the Wiedemann-Franz lawi Moreover, the 
Seebeck coefficient becomes decreased as the charge con- 
ductance increases due to the Mott relation^ Typically 
observed values of ZT are smaller than one, ZT < 1, 
and this is why thermoelectric materials have not been 
widely used in commercial applications. 

It has turned out recently, that this drawback of bulk 
materials can be overcome in nanoscale systems. Hick et 
al£ have predicted an increase of the figure of merit as 
the dimension of the system is reduced. Now, it is well 
known that the thermoelectric properties of nanoscale 
systems are strongly affected by quantum confinement 
(level quantization) and Coulomb blockade effects^— 
These effects can lead to violation of the Wiedemann- 
Franz law and failure of the Mott relation j 11 ' 12 Moreover, 
the thermal conductance of low dimensional systems is 
rather small, which allows to reach rather high values 
of ZTw^ Owing to these properties of nanoscale struc- 



tures, which are promising from the application point 
of view, interest in thermoelectric phenomena revived in 
recent years. Accordingly, thermoelectric properties of 
quantum dots (QDs), quantum wires, molecules and sil- 
icon nanojunctions have been investigated both experi- 
mentally-^— and theoretically^ 9 . - — Oscillations of the 
thermopower^ 5 - and thermal conductanc e 9 ' 1 1 with exter- 
nal gate voltage have been reported in thermoelectric 
transport through Coulomb islands. Moreover, it has 
been shown that spin correlations can strongly influence 
the thermoelectric effects in the Kondo regimei^22 Re- 
cently, an increase of the thermal efficiency has been re- 
ported in a multilevel QD for temperatures much smaller 
than the intradot charging energy^ - — Furthermore, 
the influence of ferromagnetic leads (spin polarized cur- 
rent) on thermoelectric transport properties in a single- 
level QD has been also investigated in different coupling 
rcgimesi 20 ' 21 In turn, Hatami et al. have analyzed the 
Peltier and Seebeck phenomena in magnetic multilay- 
ered structures and shown that the thermoelectric ef- 
fects significantly depend on the relative alignment of 
magnetizations in neighboring magnetic layers4i For a 
more comprehensive description of various thermoelec- 
tric transport properties in nanostructures we refer to 
recent review articles42r— 

However, the thermoelectric phenomena in multiple 
quantum dot systems are rather unexplored Espe- 
cially, double quantum dot systems may exhibit some 
new thermoelectric phenomena as they reveal a vari- 
ety of different interference effects, including Fano and 
Dicke resonance o 49 ' 50 or Aharonov-Bohm oscillations. 51 
The Fano effect appears due to quantum interference of 
waves resonantly transmitted through a discrete level and 
those transmitted nonresonantly through continuum of 
states. The effect was first observed as an asymmetric 
line shape of atomic emission spectra, but it also ap- 
pears in electronic transport through QD systems. In 
the case of double QDs, the Fano effect originates from 
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the quantum interference of electron waves transmitted 
coherently through the weakly coupled state to the leads 
and those transmitted through the state which is strongly 
coupled to the external electrodesj 49 ' 52 Very recently Liu 
et at have investigated thermoelectric effects in parallel 
double QDs attached to two metallic leads, and with a 
magnetic flux threading the quantum dot device^ They 
arrived to the conclusion that the figure of merit ZT can 
be enhanced in the vicinity of the Fano resonance. Sim- 
ilar conclusion also follows from a recent paper—, where 
the influence of electron interference in a two-level system 
on the maximum thermoelectric power is analyzed. 



II. THEORETICAL DESCRIPTION 
A. Model 

We consider two single-level quantum dots which are 
attached to normal/magnetic metallic leads. Since the 
inter-dot Coulomb interaction is small in comparison to 
the intra-dot one, the former will be omitted in the fol- 
lowing. The system is then described by Hamiltonian of 
the form 



5^ e ia d\ a d ia - t^2(d\ a d 2 a + /i.e.) + ^ U-, 



Recently, the spin voltage generated by a temperature 
gradient has been observed experimentally in a metal- 
lic magnet. 55 This novel phenomenon, called spin See- 
beck effect, offers a new way for generation of a pure 
spin current without accompanying charge current. The 
discovery of spin Seebeck effect has stimulated both ex- 
perimental and theoretical interest in the so-called spin 
caloritronics. Altough the spin caloritronics is as old as 
spin electronics^, it has been poorly investigated and 
has remained dormant for many years - except of a few 
experimental works on the thermoelectric properties of 
magnetic multilayers in the current-in-plane geometry^ 
More recently spin-dependent heat and charge transport 
in magnetic multilayers was studied experimentally in the 
perpendicular-to-plane geometry 58 ' 59 As concerns spin 
thermoelectric effects in quantum dot systems, these are 
rather poorly explored although interest in them has been 
growing recentlyj 21 ' 60 " — 



In this paper we consider thermoelectric effects in 
the system consisting of parallel coupled quantum dots 
embedded between metallic (normal or ferromagnetic) 
leads. In contrast to Ref.[53|, we include the Coulomb 
correlations. The latter are long range interactions 
and therefore can not be excluded in realistic systems. 
Moreover, one can expect a significant enhancement of 
the thermal efficiency in the Coulomb blockade regime. 
The basic thermoelectric characteristics are derived using 
nonequilibrium Green function method. The Coulomb 
correlations are taken in the Hartree-Fock approxima- 
tion and thus the formalism is relevant for temperatures 
higher than the specific temperature associated with the 
Kondo effect (Kondo temperature). Additionally, we 
investigate the spin Seebeck effect, too. In Section |TT] 
we present the model of a double quantum dot system 
attached to ferromagnetic/nonmagnetic leads and de- 
scribe basic thermoelectric phenomena to be analyzed. 
In Sec. IIIII we present the corresponding numerical re- 
sults. Finally, Sec. IIVI includes summary and general 
conclusions. 
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where the first term describes the left (j3 — L) and right 
{j3 = R) leads in the non-interacting quasi-particle ap- 
proximation. Here, c^ /3a (ck/3 CT ) is the creation (annihi- 
lation) operator of an electron with the wave vector k 
and spin a in the electrode j3, whereas £k^o- denotes the 
corresponding single-particle energy. The second term of 
Hamiltonian (JIJ describes spin conserving electron tun- 
neling between the leads and ith dot (i = 1,2), with 
Vj^ being the relevant matrix elements. The last three 
terms describe two quantum dots, where £j CT denotes the 
discrete energy level of the z-th dot, rii a — d\ a di a is the 
corresponding particle number operator, t is the inter-dot 
hopping parameter (assumed real and independent of the 
electron spin orientation), whereas E/j is the Coulomb en- 
ergy corresponding to double occupation of the i-th dot 
(< = 1,2). 

The dot-leads coupling is described usually by the 
width functions Y^- a = 2tt p%Vfvf* , where p% is the 

density of states in the lead fi for spin u. T^ a de- 
scribes the spin-dependent hybridization of the dots' lev- 
els = 1,2) and states of the /3th lead. For the sake 
of simplicity we assume that the couplings are constant 

const for 



within the electron band; T^- a {e) — Y lJO 

e G (—D,D), and r^ cr (e) = otherwise, where 2D de- 
notes the electron band width. 

The coupling parameters can be written in a matrix 
form as 



i/3 
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for f3 = L, R. To take into account various interference 
effects leading to suppression of the nondiagonal terms, 
we have introduced the parameters (?l and q^. These 
parameters are generally complex, and < \qp\ < 1. The 
spin-dependent coupling strengths can be then written 

as; rf la = r L (i ± PL ), rf 2CT = v\ la = qL r L ^(i± PL ), 

^22a = aC L (l±p L ), rf lff = aT R (l±6 PR ), rf 2(T = 
lfi CT = QrTrV^^ ± 6 PR ), and Tg a = T R (1 ± 5p R ) with 
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8 = 1 for parallel magnetic configuration and S = — 1 for 
the antiparallel one. The upper sign refers here to a 
while the lower sign refers to a =\.. Apart from this, pp 
((3 = L,R) is the polarization factor of the /3-th lead, 
Tz, and Tr are the coupling constants, and a takes into 
account difference in the coupling of a given electrode to 
the two dots. For a = the system becomes reduced to 
two QDs connected in series. 



B. Thermoelectric phenomena 

Consider first the situation when chemical potentials 
of the leads are independent of spin orientation (no spin 
accumulation). The charge current J, spin current J s , 
and heat current of electronic origin J®, flowing from 
the left lead to the right one, can be determined from 
the following formula: 



J 



-y 



h 



de 



°m J \f L (e) - fR(e)]T„{e), 
£ - fJ-L 

(3) 

where fp(s) = {exp[(e — fip)/kBTp] + 1} 1 is the Fermi- 
Dirac distribution function for the lead f3 with fip and Tjg 
denoting the corresponding chemical potential and tem- 
perature, and ks standing for the Boltzmann constant. 
Furthermore, T a (e) is the sum of transmission coefficients 
through the two conducting channels associated with two 
dots for carriers with spin a (a =t>4-)> an< 4 a = 1 for 
fj =j" and a = — 1 for a =|. T a (e) can be expressed by 
the Fourier transforms of retarded (G£) and advanced 
(G°) Green functions of the dots and by the coupling 
matrices {(3 = L, R); T a (e) = Tr[G a a T^G r a T^]. The 
Green functions have been calculated by the equation of 
motion technique in the Hartree-Fock approximation for 
the Coulomb term42 

In the linear response regime, Eq.Q can be trans- 
formed to the following formulas for charge, spin, and 
heat currents: 



(4) 
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where 8T is the difference in temperatures of the leads, 
and Sfi = eSV with SV being the voltage drop between 
the two electrodes. We remind, that according to our 
assumption, the chemical potential and temperature of 
the left electrode are /x + S/j, and T + ST, respectively, 
whereas of the right electrode are /x and T. Note, that 
in the linear response regime both Sfi and 5T tend to 



zero. In Eqs. © to $ 
integrals of the form 



L na (n = 0, 1, 2; a =t, |) are the 



de(e-fi) n ^T a (e). 
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The thermopower S is defined as the ratio of the volt- 
age drop 5V generated by the temperature difference 
6T, S — SV/5T, taken in the absence of charge current, 
J = 0. Thus, taking into account Eq.(|l]) and Eq.©, one 
obtains the following formula for the Seebeck coefficient: 
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Similarly, the charge and spin conductances, G and G s , 
can be expressed in terms of the integrals ([7]) as^ 



G = 



^ 2 E 



J 0cr, 



G s = — ^2&L 0a , 



while the thermal conductance can be written as 
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Note that the thermal conductance is determined on the 
condition of vanishing charge current. Thus, to deter- 
mine the basic thermoelectric parameters we need to find 
all the relevant integrals, Lq^j , Li a , and £20-. 

Consider now a more general situation, when spin ac- 
cumulation in the external leads becomes relevant, eg. 
due to long spin relaxation time or when an external spin 
dependent bias is applied to the system. In such a case 
we have to take into account spin splitting of the chemi- 
cal potential in the leads. In a general case, temperature 
may also be different in different spin channels. However, 
we neglect this assuming T independent of a. This may 
be justified as the energy relaxation time is much shorter 
than the spin relaxation one. Equations ([3]) take now the 
form 
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where the spin dependence of the Fermi distribution func- 
tion is now indicated explicitly. Thus, charge, spin and 
heat currents can still be expressed in the form of Eqs 
(d]) to ©, but with Sfi being explicitly dependent on a, 
Sfi — > 5 fiu, and L n<7 including now derivative of spin de- 
pendent Fermi distribution function. 

Since the bias is now spin dependent, the difference in 
chemical potentials, 6fi a , in the spin channel a can be 
written as 



8fi a = eSV a = e(SV + &SV S ), 



(13) 
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where SV is the voltage bias and 5V S is the spin bias^i. 
Of curse, SV S = in the absence of spin accumulation. 

Charge and spin currents can be written as J = GSV + 
(2e/h)G s 5V s , where G and G s are the linear charge and 
spin conductances defined above (see Eqs (9) and (10)). 
In turn, the spin current can be written as J s = G S SV + 
(H/2e)G6V s . 

The thermopower can be now calculated on the con- 
dition of vanishing simultaneously both spin current and 
charge current, or equivalent ly on the condition of van- 
ishing charge current in each spin channel. As a result, 
one can define spin dependent thermopower as 

* - w - <«> 

Equivalently, one can define spin thermopower S s 

in addition to the usual thermopower 

In turn, the heat conductance is then given by 

«=i>4i:( i2 „-£). ,i7, 

In the following we use the above formulas to calculate 
numerically the relevant thermoelectric coefficients. 



III. NUMERICAL RESULTS 

For numerical calculations we assume spin degener- 
ate and equal dot levels, £i„ = £o (f° r t = 1,2 and 
(7 =t,|)- We also assume similar magnetic electrodes, 
PL = PR = P, and symmetrical couplings, Tl = Tr = T. 
Typical experimental values of the dot-lead coupling may 
vary from a few microelectronovolts to a few milielec- 
tronovolts.— In the following it is convenient to relate en- 
ergy quantities to the dot-lead coupling. To have unique 
energy unites, we write the parameter T as T = 71^0, 
where Tq is constant and will be treated as the energy 
unit (the energy quantities will be related to T ). In turn, 
7 is a dimensionless parameter that describes strength of 
the dot-lead coupling in terms of Tq. 

The case of nonmagnetic leads is treated as a special 
case corresponding to p = 0. The parameters qh and qji 
are assumed to be real and equal, = Qr. = q- In turn, 
for the parameter a we assume a — 0.15, which indicates 
that there is a relatively large difference in the coupling 
of a given electrode to the two dots. This parameter, 
however, can be tuned by external gate voltages. For 
the sake of simplicity we also assume the same Coulomb 
parameters for the two dots, XJ\ = U2 = U. 



In a general case, the electrochemical potentials of the 
leads may be spin dependent. In other words, in ad- 
dition to the usual voltage bias one may also consider 
spin bias. The latter may appear when spin relaxation 
in the leads is slow so a spin accumulation may appear 
due to the spin bottle-neck effect (or simply when a spin 
bias is applied intentionally). This leads to a number of 
spin thermoelectric effects. We start from the simpler 
case when the spin relaxation in the leads is sufficiently 
fast to neglect spin accumulation, so the leads' electro- 
chemical potentials are independent of the electron spin 
orientation. 



A. Absence of spin bias 

Since the thermoelectric effects considered in the 
manuscript depend on spin polarization of the leads, we 
consider first the situation with nonmagnetic electrodes, 
and then proceed to magnetic leads. 

1. Nonmagnetic leads: p = 

The effects considered in this paper depend signif- 
icantly on the ratio of thermal energy and coupling 
strength of the dots to electrodes. Let us consider first 
the low temperature regime, fc^T -C T, or in our nota- 
tion, UbT/Tq <c 7. We have calculated the basic thermo- 
electric characteristics, like thermopower S, heat conduc- 
tance k, charge conductance G, and figure of merit ZT. 
Although we consider the regime fc^T <C I\ the tem- 
perature assumed here is higher than the corresponding 
Kondo temperature, so we do not take into account the 
Kondo correlations. In this transport regime one finds 
generally ZT <C 1.— However, this can be changed by 
quantum interference effects. 

Let us begin with the case of vanishing Coulomb in- 
teraction, [7 = 0. In Fig0]the thermoelectric quantities 
are plotted as a function of the dots' energy levels £0 for 
t/To = 0.8 and for indicated values of the parameter q. 
The latter parameter effectively describes strength of the 
indirect (via the leads) coupling of the dots. Such a cou- 
pling contributes (like direct inter-dot hoping term) to 
the formation of bonding and antibonding states. As a 
result, a broad peak corresponding to the bonding state 
and a narrow one corresponding to the antibonding state 
emerge in the density of states42 For q close to 1, the 
conductance peak associated with the antibonding state 
reveals the antiresonance character with a characteris- 
tic Fano line shape, whereas the conductance peak as- 
sociated with the bonding state is relatively broad and 
roughly Lorentzian. For smaller (but nonzero) values of 
the parameter q, the antiresonance is suppressed and one 
observes two peaks of different widths. Finally, for q = 
two peaks of equal widths emerge in the linear conduc- 
tance. Similar line shape and q dependence is observed 
in the electronic contribution to the heat conductance, 
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FIG. 1: (color online) Thermoelectric coefficients: (a) ther- 
mopower, (b) the figure of merit, (c) thermal conductance, 
calculated as a function of the dots' levels energy for indi- 
cated values of the parameter q. The other parameters are 
k B T/T = 0.01, t/To = 0.8, a = 0.15, U = 0, 7 = 1, and 
p — 0. The insets in (a) and (b) show respectively the ther- 
mopower S and figure of merit ZT for q = 1. 



as shown in FigfTJc). Positions of the peaks in thermal 
conductance correspond well to those in the electric con- 
ductance, similarly as in the case of a single dot^Si 

The thcrmopower S, shown in FigQJa), changes sign 
when £0 corresponds to one of the relevant resonances. 
This is a consequence of the compensation of charge cur- 
rent due to electrons by that due to holes. As a re- 



sult, there is no net charge current and no voltage drop, 
and consequently the thermopower vanishes. When the 
energy level is located below the resonance, the ther- 
mopower is negative because the majority carriers are 
holes. In turn, when the energy level is above the reso- 
nance the main carriers are electrons and thus the ther- 
mopower is positive. Note that S in Fig. 1(a) is measured 
in the units of ks /e, with e denoting electron charge 
(e < 0). However, when q is close to 1, one finds two more 
points where the thermopower changes sign. One of them 
is situated near the Fano resonance, where the conduc- 
tance vanishes due to destructive interference, whereas 
the other one is located in the valley between the two 
resonances. For other values of q only the point in the 
valley between the resonances is present, as no complete 
destructive interference occurs for smaller q. The largest 
increase of the thermopower appears in the vicinity of the 
antibonding resonance and for q = 1. This shows that 
the quantum interference has a huge impact on the ther- 
moelectric phenomena. This is also clearly visible in the 
figure of merit ZT, shown in FigQJb), which is consider- 
ably enhanced in the vicinity of the Fano antiresonance, 
where ZT is close to 1. Outside this region ZT is signif- 
icantly suppressed, even near the bonding energy level. 
When q decreases, the figure of merit ZT (which reaches 
almost unity for q = 1 near the Fano resonance) dimin- 
ishes as well. Thus, the interference effects play a crucial 
role in thermoelectric efficiency of DQD systems consid- 
ered here. However, in this coupling regime (ksT <C T) 
the Wiedemann-Franz law is not violated much as the 
Lorenz number (not shown) remains close to unity. 

Consider now the impact of Coulomb repulsion on the 
thermoelectric properties in the presence of Fano reso- 
nance, see FigO First, the Coulomb interaction leads to 
splitting of the double-peak structure (present for (7 = 0) 
in both the charge and thermal conductance, and charac- 
teristic Coulomb gaps occur— (see Fig^c) for the ther- 
mal conductance). This doubling of the resonances also 
leads to a richer structure of the thermopower. For a 
finite U the thermopower changes sign ten times, see 
FiglUa). Four points of zero S are associated with the 
four resonances located roughly at ±t and ±i — U. Next 
two such points are located near the Fano peak and its 
Coulomb counterpart, where the conductance disappears 
(so does the thermopower), whereas two other points are 
situated in the valleys between the narrow and broad 
maxima of the conductance. In the latter case the ther- 
mopower vanishes due to weighted-symmetry in location 
of the bonding and antibonding states with respect to 
the Fermi level. The current due to electrons tuneling 
through the bonding state is compensated then by the 
current due to holes tunneling through the antibonding 
level. This is a 'local' bipolar effect. The thermopower 
disappears also in the symmetry point e = —(7/2, as it 
has been explained in Ref.|2l]. However, there is one 
more point where the thermopower changes sign, namely 
this happens at the energy where a small maximum ap- 
pears in the Coulomb gap (this feature is due to a rem- 
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FIG. 2: (color online) Thermopower S (a), figure of merit ZT 
(b), and the thermal conductance (c), calculated as a function 
of the dots' levels energy for U /To = 2 and q — 1. The other 
parameters as in Fig[T] 



riant of the conductance maximum corresponding to the 
antibonding state). Correspondingly, more peaks occur 
in the figure of merit, see FigJ5Jb). From Figs 1 and 2 
also follows that the Coulomb interactions can increase 
the magnitude of ZT, which now exceeds 1. However, 
ZT is considerably enhanced only near the Fano antires- 
onances (as in the case of U = 0) and also in the middle 
of the Coulomb gap. The latter is due to a global bipolar 
effect. 

Above we have analyzed the range of ksT/T <C 1. 
The thermal efficiency, however, strongly depends on the 
ratio of thermal energy and coupling strength, and for 
ksT/T > 1 one might expect an increase of the thermal 
efficiency.^ On the other side, however, higher temper- 
ature suppresses the interference effects responsible for 
the Fano antiresonance. In Fig. [3] we show the tempera- 
ture dependence of the thermoelectric coefficients for the 
level position near the antibonding state and vanishing 
Coulomb interactions. The results clearly show that the 
thermoelectric effects are optimized for UbT/T ss 10, or 
equivalently hsT/To w IO7. To show this more explic- 




FIG. 3: (color online) Thermoelectric coefficients: (a) ther- 
mopower, (b) thermal conductance, (c) figure of merit, calcu- 
lated as a function of temperature for the position of the dots' 
levels in the vicinity of the Fano resonance, eo/Fo = —1.15, 
and for q = 1, U = 0, and 7 = 0.01. The other parameters as 
in Fig[U 



itly, in Fig. [4] we present the temperature and level po- 
sition dependence of the thermoelectric coefficients. Fig- 
ure Ufc) clearly shows that the largest enhancement of 
the thermoelectric efficiency occurs in the vicinity of the 
antibonding level. Moreover, as the temperature grows, 
a strong increase of the thermal conductance is observed 
in the valley between the two low-temperature maxima of 
the conductance, see Fig.[4|a). Such a peak is absent in 
the electric conductance and also does not appear in ther- 
mal conductance in the low temperature regime. When 
the temperature is relatively high, the Fermi-Dirac dis- 
tribution around the Fermi level becomes smeared, and 
electrons of higher energies are involved in transport. As 
the energy of tunneling electrons does not influence the 
electric conductance, it plays a crucial role in the ther- 
mal conductance. As already mentioned above, in the 
symmetry point, £0 = ~U/2, the charge current due to 
electrons is compensated by current due to holes. How- 
ever, both electrons and holes flow in the same direction, 
so the energy carried by both types of carriers is not 
compensated, in contrast to the charge. As a result, an 
additional peak centered at Eq = —U/2 appears in the 
thermal conductance for a sufficiently high temperature. 
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FIG. 4: Thermal conductance (a), thermopower (b), and the 
figure of merit (c), as a function of temperature and dots' 
levels energy for q = 1, U = 0, 7 = 0.01, and the other 
parameters as in Fig[T] 



The thermoelectric efficiency of the system under con- 
sideration can be further optimized by tuning asymme- 
try (parameter a) in the coupling of a given lead to the 
two dots. To show this, in Fig. [S] we present the a and 
level position dependence of the figure of merit for two 
distinct coupling regimes: ksT/T <C 1 (part (a)) and 
fesT/r > 1 (part (b)). The thermoelectric efficiency 
in these two transport regimes is optimized for different 
ranges of the parameter a. n the low temperature regime 
{ksT /T <C 1), ZT is optimized for intermediate values of 
the asymmetry parameter, roughly for a € (0.3,0.6). In 
turn for ksT /T > 1, the figure of merit achieves the high- 
est values for relatively high asymmetry, a £ (0.1, 0.2). It 
is worth noting that in the former case the figure of merit 
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FIG. 5: Figure of merit as a function of the asymmetry param- 
eter a and the dots' levels energy, calculated for ksT '/To = 
0.01 and 7 = 1 (a), and k B T/F = 0.1 and 7 = 0.01 (b). 
Apart from this, q = 1, U = 0, and the other parameters as 
in FigCD 



ZT can exceed 1 for properly chosen asymmetry parame- 
ter a despite of relatively strong dot-lead coupling, while 
in the latter case ZT can reach extremely large values 
exceeding 300. 

The strong dependence of ZT on the asymmetry pa- 
rameter a, displayed in Fig. [5] (a), is due to the inter- 
ference effects of electron waves passing through differ- 
ent paths in the system. More specifically, it originates 
from the a dependence of the position of the point where 
the electron conductance (as well as the thermal conduc- 
tance) reaches almost zero due to the destructive quan- 
tum interference. At this point (and in its neighborhood) 
the thermoelectric effects are maximized, which is re- 
flected by enhanced values of the figure of merit ZT in 
Fig. [5] (a). In turn, in the high temperature regime, Fig. [3] 
(b), the Fermi-Dirac distribution around the Fermi level 
becomes smeared and since the low and high energy elec- 
trons contribute differently to heat current, the thermal 
conductance peaks do not coincide now with those in the 
electric conductance (contrary to the low temperature 
regime where they do). Thus, behavior of the thermal 
conductance is now more complex. In the vicinity of the 
antibonding state, the thermal conductance is strongly 
suppressed achieving local minimum. This region of sup- 
pressed thermal conductance appears in the 'tail' of the 
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FIG. 6: (color online) Thermoelectric figure of merit ZT as 
a function of the dots' levels energy, calculated for indicated 
values of the leads' polarization p in the parallel magnetic 
configuration, and for q = 1, U — 0, ftsT/To = 0.01 and 
7 = 1. The other parameters as in Fig[T] 



electric conductance peak, where G is finite. Moreover, 
position of the minimum in thermal conductance depends 
on the asymmetry parameter a, while position of the 
peak in electric conductance is rather unchanged. This 
behavior leads to a large enhancement of ZT seen in 
Fig. [HJb) for 0.1 < a < 0.2. When a increases above 
0.2, the intensity of the peak in electric conductance cor- 
responding to the antibonding state becomes suppressed, 
and this leads to suppression of ZT for a > 0.2. On the 
other hand the suppression of ZT for a < 0.1 is asso- 
ciated with a strong dependence of the minimum posi- 
tion in the thermal conductance on the parameter a for 
a 6(0,0.2). 
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FIG. 7: (color online) Thermoelectric coefficients: (a) ther- 
mopower, (b) thermal conductance, (c) the figure of merit as 
a function of the dots' levels energy, calculated for parallel 
(P) and antiparallel (AP) magnetic configurations, and for 
p = 0.4, and U/To = 2. The other parameters as in Fig[6] 



Now we consider the situation when both electrodes 
are ferromagnetic, with the corresponding spin polariza- 
tion factor p. In the following only collinear, i.e. parallel 
and antiparallel, configurations will be analyzed. 

In the coupling regime UbT/T <C 1, the Seebeck co- 
efficient and figure of merit are only weakly dependent 
on the leads' polarization. As it has been mentioned 
above, ZT is enhanced only in the vicinity of the Fano 
resonance. For U = 0, one can notice the double nar- 
row peak structure in ZT near the Fano resonance, see 
Fig|51 where ZT is shown for three different values of 
the spin polarization factor p (including for comparison 
also the case of p = 0). Between the peaks ZT reaches 
zero due to vanishing thermopower (as already explained 
above). However, an additional double peak feature in 
ZT appears for sufficiently large p, with the intensities 
(widths) increasing (decreasing) with increasing polariza- 
tion. This feature is absent for small values of p (also for 
p = 0, see the inset to Fig. 1(b)). The appearance of the 
second double peak structure can be accounted for as fol- 



lows: Since the coupling between the dots and external 
(magnetic) leads is spin-dependent for p ^ 0, the level 
widths of the bonding and antibonding states are spin 
dependent, too. This results in narrowing (broadening) 
of the conductance peak corresponding to the spin down 
(spin up) carriers. Consequently, the characteristic zero- 
conductance Fano point is achieved first in the spin-down 
channel when going towards larger negative values of £o- 
This, in turn, results in well resolved sharp features in 
the Seebeck coefficient - one coming from spin-up con- 
tribution and another one (smaller) from the spin-down 
contribution. As a result, additional double-peak struc- 
ture appears in the thermal conductivity (not shown) and 
also in the figure of merit, as in Fig|51 The role of finite 
U is similar to that described for p = 0. 

The thermoelectric coefficients depend now on mag- 
netic configuration of the system. In Fig{7] the ther- 
mopower, thermal conductance and figure of merit ZT 
are displayed for the parallel and antiparallel magnetic 
configurations and for a finite U. It is worth noting that 
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FIG. 9: (color online) Lorentz ratio as a function of the dots' 
levels energy calculated for indicated values of the leads' po- 
larization p. Part (a) corresponds to parallel configuration 
and U = 0, while part (b) to U /To — 2 and for parallel (P) 
and antiparallel (AP) magnetic configurations. The other pa- 
rameter as in Fig|S] 



FIG. 8: (color online) Thermoelectric coefficients: (a) ther- 
mopower, (b) thermal conductance, (c) the figure of merit as 
a function of the dots' levels energy, calculated for indicated 
values of the leads' polarization p in the parallel magnetic 
configuration, and for U/Tq = 0, q = 1, UbT/To — 0.1, and 
7 = 0.1 (fcsT/r = I). The other parameters as in Fig[T] 



both thermopower and figure of merit are larger in the 
antiparallel configuration, whereas the thermal conduc- 
tance is larger in the parallel configuration. The differ- 
ence between S (and also ZT) in the parallel and an- 
tiparallel configurations is clearly visible at the corre- 
sponding peaks, where it can be quite significant. In 
turn, the difference between heat conductances in both 
configurations is rather small and well resolved only in 
the vicinity of the relevant resonances. Anyway, FigJT] 
reveals the possibility of constructing a heat spin valve, 
which would be a heat analog of the usual current spin 
valve. The suppression of thermal conductance in the 
antiparallel magnetic configuration is of similar origin as 
the suppression of linear electric conductance. The later 
was accounted for by Julliere in terms of the two cur- 
rent model^, and is due to different densities of states at 
the Fermi level in both configurations for electrons of a 
given spin orientation. More specifically, in the parallel 
configuration there is one spin channel (corresponding to 



high density of states in the source and drain electrodes) 
with high electric (and also heat) conductance and one 
with low conductance (corresponding to low density of 
states in the source and drain electrodes) . In turn, in the 
antiparallel configuration both conduction channels have 
reduced conductance as now each of them corresponds to 
high density of states in one electrode and low density of 
states in the second electrode. 

Now let us discuss briefly the temperature regime 
k B T/T « 7 (k B T/T « 1)/ In Fig. [8] the thermoelec- 
tric coefficients arc displayed for indicated values of the 
leads' polarization, no Coulomb interaction, and for par- 
allel magnetic configuration. One can notice, that the 
thermal conductance increases with increasing leads' po- 
larization p. In turn, the magnitudes of the thermopower 
as well as of the figure of merit decrease with increasing 
p. Similar behavior has been also reported for a single 
quantum dot attached to ferromagnetic leads ^ 

Although the quantum interference effects become 
smeared out by increasing temperature, the largest 
changes in the thermopower for k B T /Tq ss 1 still occur in 
the vicinity of the antibonding level, where the conduc- 
tance is strongly suppressed for k B T/T <C 1. Therefore, 
the figure of merit ZT achieves there relatively high val- 
ues (above 30 in Fig. 8(c)). All this leads to violation of 
the Wiedemann- Franz law, which can be measured by the 
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FIG. 10: (color online) Thermoelectric coefficients: (a) ther- 
mopower, (b) thermal conductance, (c) the figure of merit, 
as a function of the position of the dots' energy levels cal- 
culated for parallel (P) and for antiparallel (AP) magnetic 
configurations, and for p — 0.4. The other parameters as in 
Figi] 



Lorenz ratio, L — n/GT, as shown explicitly in Fig. [9j 
where the Lorenz ratio is displayed for indicated values 
of the leads polarization p and for interacting (U ^ 0) 
and non-interacting (U = 0) cases. One can note that 
the Coulomb interactions lead to larger deviations from 
the Wiedemann-Franz law, especially in the vicinity of 
the symmetry point £o = —U/2 [Fig. HJb)]. 

Similarly as for the temperature range studied before, 
thermoelectric coefficients depend on magnetic configura- 
tion, too. By changing magnetic configuration of the sys- 
tem from parallel to antiparallel one, the thermoelectric 
efficiency increases and is comparable to that obtained for 
p = 0. This is because in the antiparallel magnetic config- 
uration both electric and thermal conductance decrease 
(see Fig. ITOT b) for the thermal conductance). However, 
the latter suffers a larger drop and thus the figure of merit 
increases. 



B. Spin bias and spin thermoelectric effects 

Assume now that spin accumulation in the external 
leads becomes relevant, which may happen when spin 
relaxation time in the leads is sufficiently long. In such 
a case we have to take into account spin splitting of the 
electro-chemical level in both leads. Accordingly, the bias 
is then spin dependent and the difference in chemical 
potentials, A/i^, in the spin channel a can be written as 
(see section 2 for details); A^ a = eV a — e{V + aV s ). Of 
course, the spin bias V s vanishes in the absence of spin 
accumulation. 

As we have already mentioned in the introduction, 
thermally induced spin voltage has been observed re- 
cently in ferromagnetic metallic slabs. 55 A consequence 
of this is a spin analog of the Seebeck effect, so-called 
spin Seebeck effect (or spin thermopower) . In the fol- 
lowing we present some numerical results for the spin 
thermoelectric effects in the system under consideration, 
especially for the spin thermopower and spin analog of 
the figure of merit. The latter quantity is defined as 
Z S T = (2e/h)G s S s2 T/K 

In Fig. Qj] we present the spin thermoelectric coeffi- 
cients in the case of UbT /T = 1. Since the basic features 
of the thermopower S are similar to those obtained and 
discussed earlier, we focus only on the spin thermopower 
S s , see Fig. 11 (a). The spin thermopower depends on the 
leads' polarization more strongly than the thermopower 
S. Note, 5 s vanishes in the limit of p = and 15**1 grows 
up with increasing p. Comparing Fig. [TT] and Fig. [5] one 
finds that the sign of the spin thermopower S s is oppo- 
site to the sign of the charge thermopower S. This is 
because the voltage drop induced in the spin minority 
channel is larger than that induced in the spin majority 
channel. Generally, the spin thermopower vanishes at 
the same level positions as the charge thermopower does. 
Similar behavior refers to the spin conductance shown in 
Fig. ITlT b). It vanishes for p — and grows with increas- 
ing p, being positive in Fig. 11(b). The spin analog of the 
figure of merit, Z S T, is shown in Fig. [TTT c). It obviously 
grows up with increasing p as a result of the increase 
of both \S S \ and G s and a rather weak dependence of 
the thermal conductance on the leads polarization. The 
corresponding figure of merit ZT is shown in Fig. Illf d). 
ZT calculated on the condition of vanishing currents in 
both spin channels depends on the magnetic polarization 
of the leads in a more complex way than ZT obtained 
on the condition of J = 0. Roughly speaking, it grows 
up (drops) with increasing p for £o > (eq < 0). The 
suppression of ZT corresponds to the regions where the 
thermopower S is almost independent of the spin polar- 
ization of the leads. 

Qualitatively similar behavior of the spin thermoelec- 
tric coefficients can be observed also in the low temper- 
ature regime, ksT/T <§; 1, as shown in Fig. [12] for a rel- 
atively high polarization factor, p — 0.9. The maximum 
value of spin figure of merit Z S T is now slightly higher 
than that in Fig.ll, contrary to ZT which now is signif- 
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FIG. 11: Spin thermoelectric coefficients: (a) spin thermopower, (b) spin conductance, (c) Z S T, (d) ZT, calculated as a 
function of the dots' levels energy for indicated values of the leads polarization p in the parallel magnetic configuration. The 
other parameters: k B T/T = 0.1, 7 = 0.1, t/T = 0.8, a = 0.15, U/F = 2, q = 1. 
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FIG. 12: Spin thermoelectric coefficients: (a) spin thermopower, (b) spin conductance, (c) Z S T, (d) ZT, calculated as a function 
of the position of the dot levels for indicated values of the leads polarization p in the parallel magnetic configuration. The other 
parameters: ksT/Vo = 0.01, 7=1. The other parameters as in Fig. 11. 
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icantly smaller than in the case shown in Fig.ll. Apart 
from this, one feature deserves mentioning here, i.e. the 
appearance of a negative spin conductance, which ap- 
pears in some regions of the dots' energy levels. This, in 
turn, leads to negative values of the spin figure of merit 
Z S T, see Fig. [TST c). Moreover, the spin conductance is 
now enhanced in comparison to that shown in Fig. 11(b). 

IV. SUMMARY 

We have analyzed thermoelectric effects, like ther- 
mopower, electronic contribution to heat conductance, 
and thermoelectric efficiency in a system of two cou- 
pled quantum dots. The key point of the analysis was 
the role of interference effects (especially of the Fano 
antiresonances) . We have shown that the interference 
effects can significantly enhance the thermoelectric effi- 
ciency, especially for specific thermal energy ranges. The 
thermoelectric efficiency can be additionally enhanced 
by Coulomb interactions. Thus, the thermoelectric effi- 
ciency can be controlled by temperature, Coulomb inter- 
actions and quantum interference effects. Moreover, we 
have also shown that the thermoelectric efficiency can be 
additionally controlled by asymmetry of the coupling of a 
given electrode to the two dots. Some of the parameters, 
especially coupling between the two dots and between the 
dots and leads can be tuned by external gate voltages. 
It is more difficult to control Coulomb parameter U, al- 
though it can be also tuned by changing size of the dots 
via additional gates. 

We have also considered spin thermoelectric phenom- 
ena, in particular the spin analog of the thermopower 
- the so-called spin thermopower or spin Seebeck effect. 
The latter effect occurs as a result of spin asymmetry of 
the two spin channels. From the experimental point of 
view, the spin thermopower may be observed when the 
spin relaxation time in the external leads is sufficiently 
long, so the spin accumulation may build up. We have 
also analyzed the corresponding spin figure of merit. 

As we have already mentioned above, the interference 
effects (like Fano resonance) can significantly enhance the 
thermoelectric efficiency of a device. The results pre- 
sented include only electronic contribution to the heat 
current. However, one should bear in mind that there 
is also an additional phonon term in the heat conductiv- 



ity, and this contribution to heat transport may reduce 
the thermoelectric efficiency presented above, especially 
at higher temperatures . 42 ' 66 At low temperatures, how- 
ever, this contribution is small and therefore the results 
reported in this paper present a reasonable and satisfac- 
tory description. The formula for thermoelectric figure 
of merit, including the phonon contribution, can be writ- 
ten as ZT = ZT e / (1 + Kph/n e ), where K p h is the phonon 
contribution to the heat conductance and ZT e is the fig- 
ure of merit of pure electronic origin (obtained assum- 
ing K p h — 0). Since the high thermoelectric efficiency 
corresponds to a significant suppression of the electronic 
thermal conductance near the antibonding state, it is ev- 
ident that the phonon contribution to the thermal con- 
ductance, especially at high temperatures, may play a 
significant role and can not be neglected. However, it has 
been shown recently that the phonon thermal conduc- 
tance may be remarkably reduced in properly prepared 
silicon nanowire nanojunction o 15 ' 17 ! 67 , leading to impres- 
sive values of ZT (ZT > 1) at room temperature. The 
phonon contribution to the heat conductance in quan- 
tum dot systems can also be significantly suppressed by 
appropriate design of the device 6 ^ (by creating a vac- 
uum layer). Following Ref. 6&, we have estimated the 
magnitude of ZT at ksT = 10r assuming experimen- 
tally available values of T, T ps 2.5 meV (which is ex- 
perimentally available in semiconducting quantum dots) . 
Accordingly, for T ~ 300K we find the maximum value 
of ZT to be of an order of ZT w 18. To reduce the 
phonon contribution to heat conductance one can also 
design the system using phononic crystals as the elements 
of leads 6 ^. Furthermore, one can also minimize the prob- 
ability of phonon transmission utilizing strong phonon 
reflection at the interfaces constructed from materials of 
dissimilar vibrational spectral In turn, calculating the 
phonon contribution to the heat conductance according 
to Ref. |2ll one finds ZT up to 8, even in the absence of 
a vacuum layer. 

Acknowledgments 

This work was supported by Polish Ministry of Sci- 
ence and Higher Education partly as a research project 
in years 2010 - 2013 and partly as a research project in 
years 2009-2011. 



* Electronic address: |ptroch a@amu. edu.pl 
t Electronic address: barnas@amu.edu.pl 

1 R. Franz, G. Wiedemann, Ann. Phys. (Berlin) 165, 497 
(1853). 

2 M. Cutler and N. F. Mott, Phys. Rev. 181, 1336 (1969). 

3 L. D. Hicks and M. S. Dresselhaus, Phys. Rev. B 47, 16631 
(1993). 

4 C. W. J. Beenakker and A. A. M. Staring, Phys. Rev. B 
46, 9667 (1992). 



5 Y. M. Blanter, C. Bruder, R. Fazio, and H. Schoeller, Phys. 
Rev. B 55, 4069 (1997). 

6 M. Turek and K. A. Matveev, Phys. Rev. B 65, 115332 
(2002). 

7 J. Koch, F. von Oppen, Y. Oreg, and E. Sela, Phys. Rev. 
B 70, 195107 (2004). 

8 B. Kubala and J. Konig, Phys. Rev. B 73, 195316 (2006). 

9 X. Zianni, Phys. Rev. B 75, 045344 (2007). 

X.-M. Zhang, X. Chen, W. Lu, Phys. Lett. A 372, 2816 



13 



(2008). 

11 B. Kubala, J. Konig, and J. Pekola, Phys. Rev. Lett. 100, 
066801 (2008). 

12 P. Murphy, S. Mukerjee, and J. Moore, Phys. Rev. B 78, 
161406(R) (2008). 

13 R. Venkatasubramanian, Recent Trends in Thermoelectric 
Materials Research III, Semiconductors and Semimetals 
Vol. 71 (Academic Press, New York, 2001), pp. 175201; 
G. Chen, Recent Trends in Thermoelectric Materials Re- 
search III, Semiconductors and Semimetals Vol. 71 (Aca- 
demic, New York, 2001), pp. 203259. 

14 P. Reddy, S. Y. Jang, R. A. Segalman, and A. Majumdar, 
Science 315, 1568 (2007). 

15 A. I. Hochbaum, R. Chen, R. D. Delgado, W. Liang, E. C. 
Garnett, M. Najarian, A. Majumdar, and P. Yang, Nature 
(London) 451, 163 (2008). 

16 K. Baheti, J. A. Malen, P. Doak, P. Reddy, S. Y. Jang, T. 
D. Tilley, A. Majumdar, and R. A. Segalman, Nano Lett. 
8, 715 (2008). 

17 A. I. Boukai, Y. Bunimovich, J. Tahir-Kheli, J. K. Yu, W. 
A. Goddard III, and J. R. Heath, Nature (London) 451, 
168 (2008). 

18 K. Schwab, E. A. Henriksen, J. M. Worlock, and M. L. 
Roukes, Nature (London) 404, 974 (2000). 

19 L. G. C. Rego and G. Kirczenow, Phys. Rev. Lett. 81, 232 
(1998). 

20 Y. Dubi and M. Di Ventra, Nano Lett. 9, 97 (2009). 

21 R. Swirkowicz, M. Wierzbicki, and J. Barnas, Phys. Rev. 
B 80, 195409 (2009). 

22 T. Markussen, A. P. Jauho, and M. Brandbyge, Phys. Rev. 
B 79, 035415 (2009). 

23 M. Galperin, A. Nitzan, and M. A. Ratner, Mol. Phys. 
106, 397 (2008). 

24 A. M. Lunde, K. Flensberg, and L. I. Glazman, Phys. Rev. 
Lett. 97, 256802 (2006). 

25 D. Segal, Phys. Rev. B 72, 165426 (2005). 

26 F. Pauly, J. K. Viljas, and J. C. Cuevas, Phys. Rev. B 78, 
035315 (2008). 

27 Y-S. Liu, Y. -R. Chen, and Y. -C. Chen, ACS Nano 3, 
3497 (2009). 

28 J. P. Bergfield, M. A. Solis, and C. A. Stafford, ACS Nano 
4, 5314 (2010). 

29 D. Boese and R. Fazio, Europhys. Lett. 56, 576 (2001). 

30 B. Dong and X. L. Lei, J. Phys.: Condens. Matter 14, 
11747 (2002). 

31 M. Krawiec and K. I. Wysokinski, Phys. Rev. B 73, 075307 
(2006). 

32 R. Sakano, T. Kita, and N. Kawakami, J. Phys. Soc. Jpn. 
76, 074709 (2007). 

33 T.-S. Kim and S. Hershfield, Phys. Rev. B 67, 165313 
(2003). 

34 R. Scheibner, H. Buhmann, D. Reuter, M. N. Kiselev, and 
L. W. Molenkamp, Phys. Rev. Lett. 95, 176602 (2005). 

35 R. Franco, J. Silva- Valencia, M.S. Figueira, J. Magn. 
Magn. Mater. 320, 242 (2008). 

36 M. Yoshida and L. N. Oliveira, Physica B 404, 3312 (2009). 

37 T. A. Costi and V. Zlatic, Phys. Rev. B 81, 235127 (2010). 

38 M. Wierzbicki and R. Swirkowicz, J. Phys.: Condens. Mat- 
ter 22, 185302 (2010). 

39 J. Liu, Q.-F Sun, and X. C. Xie, Phys. Rev. B 81, 245323 
(2010). 

40 D. M.-T. Kuo, Jpn. J. Appl. Phys. 48, 125005 (2009). 



M. Hatami, G. E. W. Bauer, Q. Zhang, and P. J. Kelly, 
Phys. Rev. B 79, 174426 (2009). 

G. J. Snyder and E. S. Toberer, Nature 7, 105 (2008). 
J.-S. Wang, J. Wang, and J. T. Lu, Eur. Phys. J. B 62, 
381 (2008). 

Y. Dubi and M. Di Ventra, Rev. Mod. Phys. 83, 131 
(2011). 

Z.-Y. Zhang, J. Phys.: Condens. Matter 19, 086214 (2007). 

F. Chi, J. Zheng, X. -D. Lu, K. -C. Zhang, Phys. Lett. A 
375, 1352 (2011). 

M.Wierzbicki and R. Swirkowicz, Phys. Rev. B 84, 075410 
(2011). 

G. Gomez-Silva, O. Avalos-Ovando, M. L. Ladron de Gue- 
vara, and P. A. Orellana, e-print arXiv:1108.4460v2. 

P. Trocha, J. Barnas, Phys. Rev. B 76, 165432 (2007). 
P. Trocha and J. Barnas, J. Nanosci. Nanotechnol. 10, 
2489 (2010). 

D. Loss and E. V. Sukhorukov, Phys. Rev. Lett. 84, 1035 
(2000). 

M. L. Ladron de Guevara, F. Claro, P. A. Orellana, Phys. 
Rev. B 67, 195335 (2003). 

Y. -S. Liu and X. F. Yang, J. Appl. Phys. 108, 023710 
(2010). 

O. Karlstrom, H. Linke, G. Karlstrom, A. Wacker, Phys. 
Rev. B 84, 113415 (2011). 

K. Uchida, S. Takahashi, K. Harii, J. Ieda, W. Koshibae, 
K. Ando, S. Maekawa, and E. Saitoh, Nature (London) 
455, 778 (2008); K. Uchida, T. Ota, K. Harii, S. Taka- 
hashi, S. Maekawaa, Y. Fujikawa, E. Saitoh, Solid State 
Commun. 150, 524 (2010). 

M. Johnson, R.H. Silsbee, Phys. Rev. B 35, 4959 (1987). 
J. Shi, K. Pettit, E. Kita, S.S.P. Parkin, R. Nakatani, M.B. 
Salamon, Phys. Rev. B 54, 15273 (1996). 
L. Gravier, S. Serrano-Guisan, F. Reuse, and J.-P. Anser- 
met, Phys. Rev. B 73 024419 (2006). 

L. Gravier, S. Serrano-Guisan, F. Reuse, and J.-P. Anser- 
met, Phys. Rev. B 73 052410 (2006). 

Y. Dubi and M. Di Ventra, Phys. Rev B 79, 081302(R) 
(2009). 

Y. -S. Liu, F. Chi, X. -F. Yang, and J. -F. Feng, J. Appl. 
Phys. 109 053712 (2011). 

T. Rejec, R. Zitko, J. Mravlje, and A. Ramsak, e-print 
larXiv:1109.3066"V l. 

J. M. Ziman, Electrons and Phonons. The Theory of Trans- 
port Phenomena in Solids (Oxford University Press, 1960). 
W. G. van der Wiel, S. D. Franceschi, J. M. Elgerman, S. 
Tarucha, and L. P. Kouwenhoven, Rev. Mod. Phys. 75, 1 
(2002). 

M. Julliere, Phys. Lett. A 54, 225 (1975). 

N. Stojanovic, D. H. S. Maithripala, J. M. Berg, and M. 

Holtz, Phys. Rev. B 82, 075418 (2010). 

T. Markussen, A. P. Jauho, and M. Brandbyge, Phys. Rev. 

Lett. 103, 055502 (2009). 

David M.-T. Kuo and Y.-C. Chang, Phys. Rev. B 81, 
205321 (2010). 

A. N. Cleland, D. R. Schmidt, and C. S. Yung, Phys. Rev. 
B 64, 172301 (2001). 

W. Kim, R. Wang, and A. Majumdar, Nano Today 2, 40 
(2007). 

M. Tsaousidou and G. P. Triberis, J. Phys.: Condens. Mat- 
ter 22, 355304 (2010). 



